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A hypermassive neutron star (HMNS) is a possible transient formed after the merger of a neutron 
star binary. In the latest magnetohydrodynamic simulations in full general relativity, we find that a 
magnetized HMNS undergoes 'delayed' collapse to a rotating black hole (BH) as a result of angular 
momentum transport via magnetic braking and the magnetorotational instability. The outcome 
is a BH surrounded by a massive, hot torus with a coUimated magnetic field. The torus accretes 
onto the BH at a quasi-steady accretion rate ~ IOMq/s; the lifetime of the torus is ~ 10 ms. 
The torus has a temperature J; 10^'^ K, leading to copious {uv) thermal radiation. Therefore, the 
collapse of an HMNS is a promising scenario for generating short-duration gamma-ray bursts and 
an accompanying burst of gravitational waves and neutrinos. 
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Gamma-ray bursts (GRBs) are transient astrophysical 
phenomena that emit large amounts of energy (typically 
10^^ ergs) in the gamma ray band 0. The typical time 
variability is shorter than 10 ms and the duration tjur is 
~ 10 ms-1000 s. These facts suggest that the central en- 
gine of GRBs is a stellar-mass compact object, and that 
the huge energy is supplied by converting gravitational 
binding energy into radiation. The popular theoretical 
candidate for the central engine is a rotating stellar-mass 
black hole (BH) surrounded by a massive, hot accretion 
torus (see ,1] and references therein). 

Recent observations indicate that there are at least two 
classes of GRBs: short-hard GRBs (hereafter SGRBs) 
with tduT ~ 10 ms-2 s and long-soft GRBs with t^ur ~ 2~ 
1000 s. For some long GRBs, supernovae in spiral galax- 
ies have been observed coincidently indicating that 
the central engine (stellar-mass BH plus torus) for long 
GRBs is produced through stellar core collapse of mas- 
sive stars in the star forming region of spiral galaxies. 
By contrast, associations between SGRBs and elliptical 
galaxies have been reported Since elliptical galaxies 
have not produced massive stars in the past ~ 10^" yrs, 
SGRBs are most probably not related to supernova stel- 
lar core collapse. In addition, recent observations of the 
afterglow of the SGRB 050709 rule out the presence of 
a supernova light curve and point to a binary compact 
object merger as the most likely central engine 0. 

The merger of binary neutron stars (BNSs) has been 
proposed P, Q as a candidate for SGRBs. According 
to this scenario, after the merger, a stellar-mass BH is 
formed with an ambient accretion torus of mass ~ 1- 
10% of the total. The latest general relativistic hydro- 
dynamic (GRHD) simulations (with no magnetic fields) 
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have shown that just after the merger of a BNS, cither 
a BH or a neutron star is formed [3, ^3 ■ A BH forms 
promptly if the total mass of the system, M, is larger 
than a critical mass Mthr- For mergers of nearly equal- 
mass BNSs (the most likely case according to the data of 
observed binary pulsars [Til), f^'" ^^^^ than 1% of the mat- 
ter remains outside the horizon, which is unfavorable for 
GRBs. On the other hand, for M < Mthr, a hypermas- 
sive neutron star (HMNS) forms. Here, the mass is larger 
than the maximum allowed mass for rigidly rotating neu- 
tron stars with an identical equation of state (EOS). An 
HMNS is supported against collapse mainly by rapid and 
differential rotation 12J. The temperature of an HMNS 
is high {T ^ 10^^ K) because of the heat generated by 
shocks during the merger process. Although T is high 
enough to produce a large amount of neutrinos, neutrino- 
antineutrino (j^i>) pair annihilation is unlikely to generate 
a GRB fireball (consisting of relativistic e+e" pairs and 
photons). This is because pD annihilation occurs primar- 
ily inside the HMNSs, and the available energy is thus 
transferred to baryons (see e.g., 0|). However, HMNSs 
are transient objects and eventually collapse to BHs, al- 
leviating this 'baryon loading problem.' 

The value of Afthr depends crucially on the neutron 
star EOS. Recent pulsar timing observations indicate 0| 
the existence of a neutron star of mass 2.1 ± O.2M0 (one 
a error). This measurement implies that the maximum 
mass of spherical neutron stars, Mgph, is larger than ^ 
2Mq and that stiff EOSs are favored. The latest GRHD 
simulations with stiff EOSs like the one derived in [T5| . 
in which Mgph « 2-2. 2Mq, indicate that Mthr is ~ 2.7- 
2.9Mq. Thus, an HMNS is hkely to be formed after 
a merger of BNSs of canonical mass « 2.6-2.8M0 [13 
rather than a prompt collapse to a BH. 

The simulations in 10] also show that the HMNS rem- 
nants are rapidly and differentially rotating and have 
triaxial shapes. These HMNSs are secularly unstable 
since magnetic fields, viscosity, and/or gravitational ra- 
diation will transport and/or dissipate angular momen- 
tum and may trigger gravitational collapse. Recent nu- 
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FIG. 1: Upper four panels: Snapshots of the density contours for p (solid curves) and velocity vectors. The contours are drawn 
for p = 10^^ g/cm'^ X IQ"" *' g/cm^ {i = 0-9). In the last panel, a curve with p — 10^^ g/cm^ is also drawn. The (red) circle 
near the center in last two panels denotes an apparent horizon. The scale of the velocity is indicated in the upper left corner. 
The lower four panels denote the magnetic field (contours of the toroidal component of the vector potential A^p) at the same 
times as the upper panels. The solid contour curves are drawn for — 0.8(1 — 0.1j)j4,p,max,o {i = 0-9) and the dotted curves 
are for A^p = 0.08(1 — 0.2i)A^_niax,o {i = 1-4). Here, A^p = j4^_niax,o is the maximum value of A^p at t = 0. 



merical simulations [lOj suggest that gravitational wave 
emission may trigger a collapse in ~ 50-100 ms for 
M ^ 0.9Mthr ~ 2.4-2. 6Mq (this time scale wih be longer 
for smaller M). In this case, the outcome will be a BH 
with a small disk (^ O.OIM©), which is not a good can- 
didate for the central engine of SGRBs. 

The other mechanisms which transport angular mo- 
mentum are magnetic braking Il2l [l6l a nd the mag- 
netorotational instability (MRI) [iTl Il8l|. These are 
likely to play a crucial role when the magnetic fields in 
the HMNS are large enough. Magnetic braking trans- 
ports angular momentum on the Alfven time scale [T^ 
111, TA ~ R/vA 102(B/10i5 G)-i (i?/15 km)-i/2 
(M/SMq)!/^ ms, where R is the radius of the HMNS. 
MRI occurs wherever 9^^! < 18], where is the an- 
gular velocity and w is the cylindrical radius. This in- 
stability grows exponentially with an e-folding time of 
tmri = 4{dfl/dlnw)~^ 0, independent of the field 
strength. For the HMNS model considered in this paper, 
we find tmri ~ 1 ms. When the MRI saturates, tur- 
bulence consisting of small-scale eddies often develops, 
leading to angular momentum transport on a timescale 
likely to be much longer than tmri [ig- 

To study the effect of magnetic fields, we have 
performed general relativistic magnetohydrodynamic 
(GRMHD) simulations for differentially rotating 
HMNSs 6] using two new GRMHD codes Here, 
we explore HMNS collapse further by performing a simu- 



lation with the following hybrid EOS: P = Pcoid = Kip^^ 
for p < pnuc and Pcoid = K2p^'^ for p > pnuc- Here, 
P and p are the pressure and rest-mass density. We 
set Fi = 1.3, F2 = 2.75, Ki = 5.16 x 10^'^ cgs, 
K2 = Kipli-^^ and Pnuc = 1.8 x lO^^ g/cm^. 
With this EOS, the maximum gravitational mass, 
M (rest mass, Mb) is 2.OIM0 (2.32Mq) for spheri- 
cal neutron stars and 2.27M0 (2.6OM0) for rigidly 
rotating neutron stars. These are similar values to 
those in realistic stiff EOSs We construct a 

differentially rotating HMNS with the following char- 
acteristics: M ~ 2.65M0, Mb = 2.96M0, maximum 
density pmax = 9.0 x 10"'^'' g/cm'^, angular momentum 
J = 0.82GM2/c, central rotation period P^ = 0.202 ms, 
ratio of polar to equatorial radius 0.3, and rotation 
period at the equatorial surface 5.4Pc- The rotation law 
is specified in the same way as in ;6;j with the differential 
rotation parameter A = 0.8. This HMNS is similar to 
that found in the BNS merger simulation of 

For the simulation, a hybrid equation of state P — 
-fcoid + (rth — ~ £coid) is used. Here, e is the spe- 

cific internal energy, and Pcoid and Ccoid denote the cold 
part of P and e _8J . The conversion efficiency of kinetic 
energy to thermal energy in shocks is determined by Fth, 
which we set to 1.3 to conservatively account for shock 
heating. A seed poloidal magnetic field is added to the 
HMNS by specifying the (^-component of the vector po- 
tential as A^p = ^fcCi7^max[(P — Pcut),0] where Pcut is 
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FIG. 2: Evolution of the surface density and averaged ther- 
mal energy per nucleon of the torus as functions of cylindri- 
cal radius at t/Pc = 33.4 (dashed curves), 39.3 (long-dashed 
curves), 44.5 (dotted-dashed curves), and 49.9 (solid curves). 



magnetic braking. 

The star collapses at t ~ 33Pc, forming a BH composed 
of ^ 85% of the total rest mass (third panel of Fig. 1). 
Material with high enough specific angular momentum 
remains outside the newly formed BH and forms an accre- 
tion torus. However, the torus is secularly unstable, since 
magnetically-induced turbulence transports angular mo- 
mentum outward. The growth of the BH by quasistation- 
ary accretion is followed by employing an excision algo- 
rithm . The accretion rate M gradually decreases and 
eventually settles down to M ~ IOMq/s. At i ~ 50Pc, 
the rest mass of the torus is ~ O.OSM©, and the total ac- 
cretion time is thus ^ 20Pc + 0.05Mo/M ~ 10 ms. Note 
also that a coUimated magnetic field has formed along 
the rotation axis (the rightmost lower panel of Fig. 1). 

To clarify the properties of the torus, we calculate the 
surface density E and the vertically averaged thermal en- 
ergy per nucleon, (u) (see Fig. 2). The local thermal en- 
ergy per nucleon is given by u = mAr£thj where the ther- 
mal part of the specific internal energy is £th = e — £coid, 
and where rriAr is the mass of a nucleon. (We assume that 
the torus is composed of free nucleons.) Thus we have 



0.04 times the maximum pressure and Ai, denotes a con- 
stant which determines the initial strength of the mag- 
netic fields. The value of Ai, is chosen so that the maxi- 
mum value of C = B'^/SttP at i = is 3.42 x lO'^. Here, 
S'^/Stt is the magnetic pressure. This implies that the 
typical magnetic field strength is ~ 5 x 10^^ G. Such a 
large value is chosen to save computational time. Simu- 
lations with 1.9 X 10~3 ^ C ^ 7.6 X 10"^ indicate that 
a scaling relation approximately holds for smaller seed 
fields from t = to BH formation via rescaling the time 
ast/tA- The simulation is performed with a uniform grid 
of size [N -|- 1, -I- 1) for cylindrical coordinates {w, z), 
which cover the region [0, L] for each direction. Here, L is 
chosen to be 5P {R w 2.75M = 10.8 km). The grid spac- 
ings are chosen as P/100, P/120, and i?/150 {N =500, 
600 and 750), and approximate convergence is confirmed. 
Outside the HMNS, we add an atmosphere with density 
10^ g/cm'^, which is necessary when employing conserva- 
tive schemes for the hydrodynamic equations. 

In Figure 1, we show snapshots of the meridional den- 
sity contours, velocity vectors, and magnetic field lines 
for selected time slices. Following an initial period of 
linear growth (t ^ ~ 13Pc), the toroidal magnetic 
field begins to transport angular momentum from the in- 
ner to the outer regions of the star (magnetic braking), 
inducing quasistationary contraction of the HMNS 6]. 
At t ~ I a, the growth of the toroidal magnetic field 
saturates 1161 . The subsequent evolution is dominated 
by MRI 18], which distorts the poloidal magnetic field 
lines and leads to the formation of turbulent eddies on 
a scale much smaller than R (see the second lower panel 
of Fig. 1). Because of the turbulence, the matter located 
near the stellar surface is blown outward. This expelled 
material, which is connected to the fiuid in the central re- 
gion, further winds up the field lines, inducing additional 
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where g and w* denote the determinant of the spacetime 
metric and the time component of the four velocity. The 
integrals are carried out along lines of lu^constant. Note 
that eth is zero a.t t — inside the HMNS and subse- 
quently grows due to shock heating. The typical thermal 
energy per nucleon is m « 94(eth/0.1c^) MeV/nucleon, or 
equivalently, T« 1.1 x 10i2(eth/0.1c2) K. 

Because of its high temperature and density, the torus 
radiates strongly in thermal neutrinos 20,'2lj. However, 
the opacity inside the torus (considering only neutrino 
absorption and scattering interactions with nucleons) is 
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so that the neutrinos are 



effectively trapped 21]. Here, T12 = K. Since 

the torus is optically thick, the neutrino luminosity may 
be estimated in the diffusion limit as L^, ^ nR'^F, 
where R is the typical radius of the emission zone, the 
flux is P - (c/3)(7iV^/4)(crrVr) (where a is the Stefan- 
Boltzmann constant and N^, is the number of neutrino 
species, taken as 3), and the neutrino optical depth is 
T - kS. Then - 2 x 10^^ erg/s(P/10 km)^T^^j:^^, 
which is comparable to the neutrino Eddington luminos- 
ity ^1]. Here, Eis = E/IO^^ g cm"^ Because of the 
geometry of the torus, pair annihilation will be most ef- 
ficient near the z-axis. Furthermore, the surface den- 
sity along the z-axis (w = 0) outside the apparent hori- 
zon is E18 ^ 0.01 for t 40Pc, which is much smaller 
than the surface density of the torus. In fact, the total 
mass contained in a cylinder of radius tu ^ AI (~ 5km) 
is ~ 1O~^M0, which is likely small enough to allow the 
formation of a relativistic fireball |1| . 
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Our numerical results suggest the presence of a hot, 
hyperaccreting torus which is optically thick to neutri- 
nos. A model for the neutrino emission in a similar flow 
environment with comparable Ly (a 'neutrino dominated 
accretion flow') is provided by Di Matteo et al. j23|. Ac- 
cording to this model, the luminosity due to vv annihi- 
lation is L^jj ~ 10*^° ergs/s Aloy et al. simulate 
the propagation of jets powered by energy input along 
the rotation axis (as would be supplied by the uv' anni- 
hilation). They find that if the half-opening angle of the 
energy injection region is moderately small 45°) and 
the baryon density around the BH is sufhciently low, jets 
with the Lorentz factors in the hundreds can be produced 
given an energy input L^y ^ 10'** ergs/s lasting ~ 100 
ms. They also show that the duration of SGRBs may 
be 10 times longer than the duration of the energy in- 
put because of the differing propagation speeds of the jet 
head and tail. Our numerical results, along with the ac- 
cretion flow and jet propagation models of (21|,|23|, thus 
suggest that magnetized HMNS collapse is a promising 
candidate for the central engine of SGRBs. Since the 
lifetime of the torus is ~ 10 ms in our simulation, the to- 
tal energy of the vv annihilation {E^^ ~ 10*^ ergs) may 
be sufficient to power SGRBs as long as the emission is 
somewhat beamed (and beaming is probably encouraged 
by the fat geometrical structure of the torus '2^). 

Alternatively, a relativistic outflow could also be pow- 
ered by MHD effects [l|. Though the GRMHD equa- 



tions are solved self-consistently in our simulation, we do 
not find evidence for strong MHD outflows. This may 
be a consequence of our initial magnetic field configu- 
ration or our neglecting neutrino pressure, and requires 
further study. However, simulations of magnetized accre- 
tion tori in fixed Kerr spacetime have found outgo- 
ing electromagnetic energy due to the Blandford-Znajek 
(BZ) effect 25]. The BZ luminosity 24J is estimated as 
Lbz 1053a2(B/10i6 G)2(M/2.8Mo)2 erg/s, where a 
is the nondimensional BH spin parameter and B is the 
typical magnetic field strength. Assuming a reasonable 
conversion efficiency from the Poynting flux to the kinetic 
energy of the fireball and then to gamma-ray energy, en- 
ergy fluxes of this magnitude are sufficient for forming 
SGRBs. 

Finally, this model predicts that SGRBs should ac- 
company a burst of gravitational radiation and neutrino 
emission from the HMNS delayed collapse. We plan to 
study this gravitational radiation in a future paper 26] . 
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